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The origin of the work function 
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In this paper we analyse the mechanisms responsible for the bonding of electrons to metal surfaces. 
We present and validate a method to measure the energy distribution of dense electron ensembles 
at ambient conditions. The method relies on the detection of the density of states near the Fermi 
level. We have found sharp structures in the energy distribution of charge confined at the surfaces 
of Zn, Cu, Al, Si and C. These structures correspond to the onset of space charge i.e. electron 
emission. The threshold for emission is at an energy level that is associated to the work function <j> 
of the metal under consideration. The work function increases linear with the free electron density 
and equals the exchange energy given by the Hartree-Fock model of electrons. We conclude that 
electron-electron interaction is the major factor that contributes to the observed work function of 
metals. In our experiments, electron emissions are triggered by the photoelectric effect. Remarkably, 
the emissions that have been observed do not obey the Einstein-Millikan threshold relation h nu=(j> 
of the photoelectric effect. By discrimination between requirements for energy and momentum we 
deduce that the emissions are trigged by photon momentum. We claim that the photoelectric effect 
is the result of a phase transition in an electron gas that connects a condensed surface state to a free 
continuum state. The assignment of a phase transition to the effect under consideration is justified 
by analysing the behaviour of the Fermi-Dirac energy distribution upon raising the Fermi level. 
When the Fermi level crosses the continuum level we predict an abrupt increment in the number of 
continuum state electrons that exceeds 70 orders of magnitude. Both the observed threshold energy 
and the line shape of the transition are in good agreement with theoretical predictions. 
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I. INTRODUCTION 

Quantum electronics has been an important research 
topic in industry and academia for the last century. Re- 
cent work includes studies on devices for cold electron 
emission using nano structures, studies on photonic de- 
vices and (bio-) chemical sensors [1-9]. Major efforts have 
been put to understand the specific role of surface states 
[9-14]. More fundamental issues concerning quantum 
many-body theory are subject to advanced experimental 
research using dense particle systems that are confined in 
optical lattices [15-17]. In the latter studies ions, atoms 
and molecules have been used to mimic degenerate elec- 
tron gases. This approach allows the full control over 
many features of interest e.g. lattice parameters or dif- 
ferences caused by the symmetry of the wave function of 
ensembles, which accounts for the statistical behaviour 
of particles. This kind of experiments has been used for 
fundamental tests on the interaction between fermions, 
bosons and mixtures. Ultimately, this will lead to a full 
understanding of a long-standing problem known in solid 
state physics: the legitimacy of the independent electron 
approximation [18]. 

Whereas the majority of the experiments have been 
carried out under vacuum conditions at low tempera- 



tures [1-17] we will study dense electron systems at am- 
bient conditions. Electron densities commonly found in 
metals yield ~ 10 23 cm~ 3 . The corresponding deBroglie 
wavelength of electrons near the Fermi level exceeds the 
mean electron separation (Table I), even at a tempera- 
ture of 300K. Thus, at ambient conditions the electron 
system in metals has to be treated as a quantum system 
in the low temperature limit [20, 21]. In this regime of 
high density and low temperature one has to go beyond 
the independent electron approximation and account for 
electron-electron interactions [18,21]. 

We consider the Hartree-Fock theory of the free elec- 
tron gas [18] as a theoretical framework for the evalua- 
tion of exchange mechanisms. Although the many elec- 
tron wave function of an ensemble is treated by a single 
particle wave function this theory accounts for electron- 
electron interactions. The electron-electron interaction 
gives rise to the so-called exchange energy term that is 
given as a function in the electron density (1/V S ) by [18] 
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This term represents the average binding energy of 
a single electron to the ensemble as a result of ex- 
change interaction. Note that the attractive coulomb 
interaction [21,22] of the proton-electron system (hy- 
drogen) is identical to the exchange interaction of the 
electron-electron system, except for an (exact) factor of 
(#t) 1/3 =0.4562720. 



2 



workfunclion data ref [18] + 
workfunclion data ret [29] x 
Hartree Fock exchange energy 
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FIG. 1: The work function (data points) and the Hartee-Fock 
estimate of exchange energy due to electron-electron interac- 
tion (solid line) for metals with different electron densities. 



Evidently, the exchange mechanism is the origin of an 
attractive coulomb force between electrons, which is re- 
sponsible for the cohesion of dense electron ensembles. 
In Fig. 1 we have plotted the exchange energy as a func- 
tion of the electron density for the metal elements. We 
have made a comparison with the work function of these 
elements known from literature [18,29]. 

From the many observations that have been made e.g. 
by photon spectroscopy and thermo ionic emission [18,19] 
it is apparent that electrons are bound to a metal by some 
kind of process. However, the origin of the bonding pro- 
cess is not well understood [18,19] . The work function has 
been originally introduced as a phcnomcnological quan- 
tity to account for the bonding of electrons to metals. It 
is generally assumed that ion-electron interactions at the 
metal-air interface play an important role in the binding 
process [18]. Based on the analysis presented in Fig. 1 
it is concluded that the work function is correlated to 
exchange energy i.e. we suggest that the binding is more 
electron-electron like. This is similar to the original work 
of N.D. Lang and W. Kohn [22] at the start of the for- 
mulation of density functional theory [18]. 

However, in this paper we will investigate this hypoth- 
esis more carefully using statistical physics and exper- 
imental proofing rather than attempting to understand 
the physics by ab inito calculations as usually explored in 
quantum chemistry. In our experiments we shall use the 
dense electron ensembles that are already present at the 
surface of common metals. We shall consider the valence 
electrons as free electrons near the Fermi level. In a next 
step, we will introduce a method to measure the density 
of states of an electron ensemble at the Fermi level. By 
observing the structure in the density of states of an en- 
semble as a function of the number of electrons present 
we shall be able to quantify the contribution of exchange 
energy to the work function. 



II. THEORETICAL DESCRIPTION 

In our experiments we will consider surface and space 
charge and treat them as free electron gases. We will con- 
trol the Fermi level of ensembles by changing the total 
number of particles that is contained in the gas. There- 
fore, we will first evaluate the statistical aspects of a lo- 
calised ensemble that contains a variable amount of elec- 
trons. The electrons present in a gas obey Fermi-Dirac 
statistics and the energy distribution is given by [18-20] 

FD(E,V)= e{E _ Ef 1 )/kT + i (2) 

for E > Ef and where Ef denotes the Fermi energy of 
the ensemble. The continuum energy level is located at 
E = 0. The total number of electrons that is present is 
determined by the external static electric potential of the 
ensemble which is denoted by V. The Fermi level is given 
by Ef — Ef — qV where Ej a is the Fermi level in ab- 
sence of this external potential. At a negative potential 
the Fermi energy is increased as a result of the addi- 
tion of charge to the system. When electrons are added, 
the Fermi level rises to the value Ef Q — qV. This is re- 
quired because the number of electrons in the ensemble 
has increased whereas the maximum occupation number 
of successive energy states is limited to two (Pauli Ex- 
clusion Principle). Thus, the Fermi level of an electron 
ensemble is shifted in presence of an external potential 
by an amount of — qV whereas the ground state level is 
fixed. 

Now, we introduce a potential well to represent a con- 
fining structure (e.g. a crystal). As usual [18,20], we 
define a constant potential of magnitude —<p — Ef in the 
region where the crystal is present and zero elsewhere. By 
introducing this potential both the ground state level and 
the Fermi level in the crystal are shifted by an amount 
—(j) — E} [18-20]. Both levels are shifted by the introduc- 
tion of a confining potential but the number of electrons 
did not change. This means that the Fermi energy has 
not changed. Thus, the Fermi level is now located at 
— <fi — qV and the ground state level at —cf)~Ef . 

By choice of these definitions, electrons are in the con- 
tinuum whenever E > 0. Secondly, the work function 
is defined by <j). The work function corresponds to the 
minimum energy that is required to promote an electron 
originally at the Fermi level of a charge neutral metal 
into the continuum [18-20]. 

The energy levels of the idealised electron system un- 
der consideration are schematically depicted in Fig. 2. 
This is the basic model that is used throughout this pa- 
per to clarify how the energy distribution of an electron 
ensemble is affected by charging a metal. In Fig. 4 we 
have evaluated the energy distribution of an electron en- 
semble held by a Zn plate (E fo = 9.5eV, (j> = A.2eV) for 
the case that the plate is grounded, and for the case that 
it is charged to — 5V. 
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TABLE I: Electronic properties of different metals relevant to this work [18]. 



Element 


Density 


Mean 


Fermi 


deBroglie 




Work 






separation 


level 


wavelength 




function (eV) 




(xl0 22 crrT 3 ) 


(r s /a ) 


(eV) 


(nm) 


[29] 


[18] 


Be 


24.7 


1.87 


14.3 


0.31 


4.98 


3.92 


Al 


18.1 


2.07 


11.7 


0.34 


4.25 


4.25 


Fe 


17 


2.12 


11.1 


0.35 


4.6 


4.31 


Mn 


16.5 


2.14 


10.9 


0.35 


4.1 


3.83 


Sb 


16.5 


2.14 


10.9 


0.35 


4.63 


4.08 


Ga 


15.4 


2.19 


10.4 


0.36 


4.2 


3.96 


Sn 


14.8 


2.22 


10.2 


0.36 


4.42 


4.38 


Bi 


14.1 


2.25 


9.9 


0.37 


4.22 


4.4 


Zn 


13.2 


2.3 


9.47 


0.38 


4.05 


4.24 


Pb 


13.2 


2.3 


9.47 


0.38 


4.25 


4 


In 


11.5 


2.41 


8.63 


0.4 


4.12 


3.8 


Tl 


10.5 


2.48 


8.15 


0.41 


3.84 


3.7 


Cd 


9.27 


2.59 


7.47 


0.43 


4.22 


4.1 


Hg 


8.65 


2.65 


7.13 


0.44 


4.49 


4.52 


Mg 


8.61 


2.66 


7.08 


0.44 


3.66 


3.64 


Cu 


8.47 


2.67 


7.0 


0.44 


4.64 


4.4 


Au 


5.9 


3.01 


5.53 


0.49 


5.31 


4.3 


Ag 


5.86 


3.02 


5.49 


0.5 


4.54 


4.3 


Nb 


5.56 


3.07 


5.32 


0.5 


4.33 


3.99 


Li 


4.7 


3.25 


4.74 


0.53 


2.9 


2.38 


Ca 


4.61 


3.27 


4.69 


0.54 


2.87 


2.8 


Sr 


3.55 


3.57 


3.93 


0.59 


2.59 


2.35 


Dr. 

Da 


o. 10 


O. / 1 




U.Ol 


9 7 


9 AQ 


Na 


2.65 


3.93 


3.24 


0.65 


2.75 


2.35 


K 


1.4 


4.86 


2.12 


0.8 


2.3 


2.22 


Rb 


1.15 


5.2 


1.85 


0.86 


2.16 


2.16 


Cs 


0.91 


5.62 


1.59 


0.92 


2.14 


1.81 


n-Si 


5.2 


3.14 










C (graphite) 


~ 3 10 18 cm -3 


81 











When a negative potential is applied the charge den 
sity has increased. For applied voltages less than —4.2V 
the Fermi level exceeds the continuum level. All of th( 
electrons that are added beyond this point are free fron 
an energetic point of view. These electrons should, ii 
principle, be released from surface into the continuum 
In our experiments we will verify that this is actually th< 
case by charging metals plates and probe for the emis 
sion of electrons. Before doing so, we have to consider th( 
density of states of surface charge and relate the change: 
to the macroscopically observable charge and voltage. 

We argue that the correct quantum mechanical de 
scription of surface charge is given by an ensemble o 
free, independent electrons with two degrees of freedom 
In classical electrostatic theory it is well established that 
negative charge at a conductor is located at the outer 
surface [23]. This is a direct consequence of Gauss' law 
that implies that in absence of an internal electric field 
the surfaces of conductors are equipotentials. The energy 
distribution of surface charge is given by the definition 
of capacitance C through Q = CV [22]. As the electric 
charge is quantized by the electron charge — q we can 
write Nq = —2CV for a macroscopic amount of charge 
that contains N electrons. A factor of two is added to 
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FIG. 2: Relevant energy levels of a degenerate electron sys- 
tem confined in a potential well. The number of electrons is 
varied by means of an external potential. In absence of such 
a potential the levels are filled over a range of energies corre- 
sponding to the Fermi energy Ef a . The Fermi level is located 
at —<f>. Electrons are added to the system by charging the 
system to a negative potential V . After charging the Fermi 
level has risen to the level — <j> — qV. 
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FIG. 3: The energy distribution of electrons at a Zn surface 
with 0=4. 2eV and Ef a =9.7eV at a bias potential of OV^ (solid 
line) and for — 5V (dashed line). Beyond —4.2V the Fermi 
level is extended into the continuum (Fermi-Dirac evapora- 
tion) . 



For an electron in the ground state (i=l, j=l), L 2 can 
be eliminated by combining Eq. [5] and [6] to obtain the 
density of states (dos) in terms of surface charge density 



dn a 



dE 



1 da 
qdE 



m 



(7) 



This result was obtained after quantisation of the free 
particle wavefunction on a square. However, even on a 
flat surface the coulomb force is a central force. Thus, 
in principle the quantisation procedure should be car- 
ried out on a disk. This involves the introduction of 
Bessel equations and its solutions. Its introduction leads 
to mathematical complexities and does not lead to more 
insight in the underlying physics. Alternatively, we adopt 
the direct counting procedure of electron states as em- 
ployed in standard textbooks on quantum physics [18-21]. 
This counting procedure is well established and involves 
a quantisation procedure that is effectively carried out 
in spherical coordinate system. The surface density of 
states that is obtained in this approximation is given by 



account for two spin states. If we use the definition of 
the Volt, which is the energy per unit of charge given by 
E = —qV, and differentiate, we obtain the energy rela- 
tion for the classical surface density of electrons which is 
given by 



dn 



dn c 

He 



2C 



(3) 



Here A denotes the surface area under consideration. 
In addition, it is well known that the capacity of metal 
objects is merely independent of the voltage if the elec- 
trical field strength just outside the object is limited to 
20 kV/cm (in order to prevent electrical breakdown in 
air). To summarize: we have identified the capacity as 
a measure of the classical surface density of electrons i.e. 
the number of electrons per unit of energy. The quantity 
is independent of energy. 

As a quantum mechanical description of surface charge 
we consider a free electron gas with two degrees of free- 
dom. In general, the single particle density is given by 
the expectation value divided by quantisation volume Q 



1 , , 2 1 



(4) 



For a free electron that is confined at a square surface 
of dimensions L x L, the charge density a is given by 



CT = JL |$| 2 = JL 
L 2 L 2 



The electron energy is given by [18-21] 



Ej 



h 2 7T 2 

2mL 2 
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(5) 



(6) 



qm 



dE 



m 



(8) 



The expectation value of the particle density of a sur- 
face quantum gas is independent of its energy similar as 
in the classical description of surface charge. Upon com- 
parison of the classical result with the quantum mechan- 
ical result we can define a minimum length scale where 
quantum mechanical effects are to be anticipated. The 
macroscopic electron surface density obtained by Gauss' 
law can be expressed as 



dn c 
~dE 



2C_ 
Aq 2 



2eo 
q 2 r 



(9) 



where we have used C = Aire r for a classical sphere 
of radius r. When we demand that this is equal to the 
density of states of the quantum mechanical electron gas 
we obtain an expression of the radius of a charged sphere 
that carries the charge of a single electron. This radius 
is given by r = 2ire h 2 jmq 2 which is equal to half of 
Bohrs' radius for hydrogen. For charge densities larger 
than 2q/a Q only a quantum mechanical description of 
the electron is appropriate. For charge densities smaller 
than 2q/a Q classical electrostatic theory is expected to 
hold. Note that this is consistent with the estimate of 
exchange energy in the Hartree-Fock model: exchange 
effects are only relevant at densities in the order of q/a 
and vanish at low density. 

Apparently, a classical plane surface relates to a topo- 
logical structure that does not contain microscopic struc- 
tures smaller than a . Thus for classical electrostatic the- 
ory to hold, a metal surface should be free of quantum 
dots (zero Degrees Of Freedom, DOF), quantum wires 
(one DOF) or quantum wells (two DOF). This is consis- 
tent with the fact that the density of states of latter struc- 
tures is dependent on energy [18,24]. Its presence leads 
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to distinct structures in the observed density of states 
(often referred to as quantum capacitance) as measured 
by e.g. conductance spectroscopy [24]. 

That the equivalency of Eq. 8 and 9 makes sense from 
a physical point of view can be demonstrated more gener- 
ally by using the virial theorem [21]. This theorem states 
that the average kinetic electron energy of a free electron 
is twice the average action of the (repulsive) coulomb 
force: 



r n 2 n 2 k 2 

2E k = / F(r')r'dr' = = — (10) 

J r >>r Aite r m 



which is implicated by setting Eq. 8 equal to Eq. 9. Note 
that the path integral of the action only contributes in the 
outer region of a charged sphere: the contribution over its 
interior is zero. From this result we deduce that a 'free' 
electron is always under influence of its own repulsive 
coulomb force [25] and that the (positive) coulomb energy 
is associated to its kinetic energy. 

The eigen energies of a quantum gas of independent 
surface electrons are given by Eq. [(51 In Table II we 
have explicitly given the lowest lying single particle en- 
ergy levels. The levels are enumerated according to the 
momentum eigen states Pij = hk(i,j) with k = ir/a and 
where a is the quantisation length of the problem under 
consideration. In Fig. 4 we have plotted the density of 
states as a function of the electron energy for a gas that 
contains 10 4 electrons. We have taken a = 1.0 10~ 7 m 
and have assumed an energy resolution of 25 meV. This 
corresponds respectively to the typical surface roughness 
of the metal surfaces and to the typical energy resolution 
in the experiments that are discussed in this paper. From 
Table II and Fig. 4 we learn that the charge density per 
unit of energy is a constant, as it should. However, this 
is only true when we consider the average of a sufficient 
number of states. 

Fluctuations originate from changes in the level of de- 
generacy of successive energy levels (see Table II). It is 
interesting to note that these fluctuations are directly 
related to the mathematical problem known in number 
theory on the representation of a natural number as the 
sum of two squares. 

Based on the discussion in the previous paragraphs 
we summarize that surface charge as treated by classical 
electrostatic theory is the limiting case of a free, quantum 
electron gas with two degrees of freedom. The classical 
limit corresponds to the quantum mechanical result in 
the limit of large quantisation lengths (or equivalent in 
the limit of H — > 0) as required by Bohrs' correspon- 
dence principle [21]. We shall confirm experimentally 
that the charge density is independent of energy and that 
the characteristic fluctuations shown in Fig. 4 are real. 
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FIG. 4: Calculation of the charge density as a function of 
energy for a two dimensional free electron gas that contains 
10,000 electrons. 



III. EXPERIMENTAL SET-UP 

We use polycrystalline plates of Zn, Cu, Al and C (as 
graphite) with dimensions of 40x40x2mm and a n-Si(100) 
plate as an electrode of an air capacitor. An aluminium 
grid (70% open, 1 mm mesh) forms the counter electrode. 
We study the transient current and voltage characteris- 
tics of this system by means of a discharge through a 
resistor. We use the photoelectric effect to liberate elec- 
trons from the surface. Emissions are observed at nearly 
field free (dV/dx<10V/cm) conditions. Experiments are 
carried out under ambient conditions. That is, at a pres- 
sure in the range of 950-1050 mbar, at temperatures in 
the range of 20-22°C and in an open-air environment. 
The electrodes are continuously flushed with nitrogen (~ 
2 L/min) to prevent oxidation. The experimental set-up 
is shown in Fig. 5. 

We use a cold cathode UV-C source (Philips UK TUV 
PL-S 9W) that emits up to 253 nm or a 150W quartz Hg 
lamp that emits up to 180 nm. The lamp spectrum of the 
150W Hg lamp is given in Fig. 6, which was measured 
using an Avaspec 2048-2 spectrometer supplied with a 
UA-178-1100 grating. Both lamps are operated by stan- 
dard lamp gear at a line frequency of 50 Hz. The light 
intensity at the surface of the metal plate is characterized 
by the light intensity in the 200-400 nm wavelength re- 
gion and was measured by a SG01S SiC UV photo diode 
(Cree Research, US). The UV intensity at the surface is 
typically 24mW/cm 2 (150W lamp) and 5.6mW/cm 2 for 
the PL-S lamp. 

The voltage (V) and current (/) as indicated were 
measured using a Keithley 6 10C electrometer at an 
impedance setting of lOOGOhm. Alternatively, the en- 
tire electro meter is replaced by a high impedance volt- 
age divider (100:1) using two resistors (Ohmcraft HVRW 
lOOGOhm and HVRW IGOhm). By doing so we have 
verified that the voltage readings are not affected by bi- 
ases that may originate from amplifier offset or unwanted 
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TABLE II: Lowest lying energy eigen states of an isotropic 

free particle quantum gas with two degrees of freedom. The 

momentum eigen states are given by pij = hk(i,j) at energy 

# i j i 2 +f #(cont) i j i 2 +f 

36 6 1 37 

1 1 1 37 2 6 40 

2 1 1 38 6 2 40 

3 1 1 2 39 4 5 41 

4 2 4 40 5 4 41 

5 2 4 41 3 6 45 

6 1 2 5 42 6 3 45 

7 2 1 5 43 7 49 

8 2 2 8 44 7 49 

9 3 9 45 1 7 50 

10 3 9 46 5 5 50 

11 1 3 10 47 7 1 50 

12 3 1 10 48 4 6 52 

13 2 3 13 49 6 4 52 

14 3 2 13 51 7 2 53 

15 4 16 52 3 7 58 

16 4 16 53 7 3 58 

17 1 4 17 54 5 6 61 

18 4 1 17 55 6 5 61 

19 3 3 18 56 8 64 

20 2 4 20 57 8 64 

21 4 2 20 58 1 8 65 

22 5 25 59 4 7 65 

23 3 4 25 60 7 4 65 

24 4 3 25 61 8 1 65 

25 5 25 62 2 8 68 

26 1 5 26 63 8 2 68 

27 5 1 26 64 6 6 72 

28 2 5 29 65 3 8 73 

29 5 2 29 66 8 3 73 

30 4 4 32 67 5 7 74 

31 3 5 34 68 7 5 74 

32 5 3 34 69 4 8 80 

33 6 36 70 8 4 80 

34 6 36 71 9 81 

35 1 6 37 72 9 81 



contact potentials. Data are logged using a 16 bits ADC 
(Meilhaus FS 1608) at a rate of 6.6 samples per second. 
The accuracy of the voltage readings is 25mV and the ab- 
solute value is accurate within ±100 mV. The potential 
of the metal plate under consideration is initially set us- 
ing an auxiliary voltage source. Slowly varying voltages 
and currents are recorded when the system is allowed to 
discharge through the impedance of the electrometer. 

Oxide layers on the metal plates are removed using 
sanding paper (P400). The final surface roughness that is 
obtained after this procedure is 0.1-1.0 \xm. By applying 
this procedure we obtain a randomly oriented polycrys- 
tallinc surface. The silicon plate is a phosphorus doped Si 
wafer (n-Si, 5 ft cm, oriented in the (100) direction, Ok- 
metic, Finland). It has an optical grade surface (surface 
roughness 1-3 nm) and is not scoured. 

In order to be able to measure the energy distribution 



uv 



10 mm 




FIG. 5: Schematic representation of the experimental set-up 
(side view). Coupons of Zn, Cu, Al, C or Si with dimensions of 
40x40x2 mm are contained in a grounded shielding box. The 
top of the box is covered by an aluminium grid. Photons are 
introduced from the outside as depicted. A battery is used 
to provide a fixed bias voltage. Initially the metal plate is 
negatively charged by an auxiliary voltage source (not shown). 
Experiments are done at ambient conditions. 
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FIG. 6: Emission spectrum of the 150W quartz mercury lamp. 



over an energy range of -42eV to 50eV we initially charge 
the metal plate to -50V to hold an excess amount of elec- 
trons. The capacitor is then discharged by the RC circuit, 
either in the dark, or under illumination with UV light. 
A battery pack of +42V is used to provide the required 
positive bias and a lOnF capacitor is used to increase the 
amount of charge in the system that is initially present 
(Fig. 5). 



IV. OBSERVATIONS AND ANALYSIS 

In absence of UV light the system discharges 
through the impedance of the voltmeter according to 
V(t)=V(0)e-*/ flc ' with R=100Gf7, C^lOnF and V(0)=- 















36 


6 


1 


1 





1 


1 


37 


2 


6 


2 


1 





1 


38 


6 


2 


3 


1 


1 


2 


39 


4 


5 


4 





2 


4 


40 


5 


4 


5 


2 





4 


41 


3 


6 


6 


1 


2 


5 


42 


6 


3 


7 


2 


1 


5 


43 





7 


8 


2 


2 


8 


44 


7 





9 





3 


9 


45 


1 


7 


10 


3 





9 


46 


5 


5 


11 


1 


3 


10 


47 


7 


1 


12 


3 


1 


10 


48 


4 


6 


13 


2 


3 


13 


49 


6 


4 


14 


3 


2 


13 


51 


7 


2 


15 





4 


16 


52 


3 


7 


16 


4 





16 


53 


7 


3 


17 


1 


4 


17 


54 


5 


6 


18 


4 


1 


17 


55 


6 


5 


19 


3 


3 


18 


56 





8 


20 


2 


4 


20 


57 


8 





21 


4 


2 


20 


58 


1 


8 


22 





5 


25 


59 


4 


7 


23 


3 


4 


25 


60 


7 


4 


24 


4 


3 


25 


61 


8 


1 


25 


5 





25 


62 


2 


8 


26 


1 


5 


26 


63 


8 


2 


27 


5 


1 


26 


64 


6 


6 


28 


2 


5 


29 


65 


3 


8 


29 


5 


2 


29 


66 


8 


3 


30 


4 


4 


32 


67 


5 


7 


31 


3 


5 


34 


68 


7 


5 


32 


5 


3 


34 


69 


4 


8 


33 





6 


36 


70 


8 


4 


34 


6 





36 


71 





9 


35 


1 


6 


37 


72 





9 



7 




-60 1 1 1 1 1 1 1 1 

200 400 600 800 1000 1200 1400 

time (s) 

FIG. 7: Measured discharge currents in presence and absence 
of the UV light. The photoelectric effect is responsible for 
deviations in the exponential decay of the discharge current. 

92V. In presence of UV light, the capacitor is additionally 
discharged by a photoelectric current. Typical discharge 
curves that have been observed are plotted in Fig. 7. 

For a quantitative study on the underlying physical 
processes we study the density of states near the Fermi 
level. We therefore analyse and interpret our data in 
terms of changes in charge (dQ) and changes in energy 
(dE) rather than by current (/) and voltage (V). The 
conversion of the variables (V,/) into (Q,dQ/dE) is ac- 
cording to E — —qV and 

dQ(E) _ 1 I(t)dt _ 1 1{t)dt 
dE ~ qV(t + dt) - V{t) ~ ~qdV{t) ' ' ^ 

This step is realised in our analysis by numerical differ- 
entiation. Note that the total charge per unit of energy 
dQ(E) / dE ~ qdn(E)/dE is proportional to the density 
of states but requires the definition of a specific surface 
area or volume (sec section V). 

Our expression for the density of states is different from 
the one that is commonly used in scanning tunneling 
spectroscopy [24]. In the latter case, the local density of 
electron states is assumed to be proportional to dl/dV. 
we believe, however, that the density of states is a static 
parameter. The use of the differential conductance as a 
direct measure of the density of states is incorrect from 
a dimensional point of view. 

In Fig. 8.1a-8.5a we present the measured density of 
states of surface charge as a function of energy for Zn, 
Cu, Al, C and n-Si in absence of light. We have observed 
a constant level on average i.e. independent on energy. 
The fluctuations in the observed charge density have been 
attributed to the fluctuations also observed in Fig. 3. In 
our experimental set-up there are two contributions to 
these fluctuations. One from the metal foil contained 
in the lOnF capacitor and a second contribution is from 
the metal-air-aluminium grid capacitor (~ lpF) under 
consideration. 



In our experiments we occasionally observe structures 
in the density of states at specific energy levels. We have 
found that these structures are the result of electronical 
surface processes (Redox reactions) at the metal-air inter- 
face. These structures appear as peaks upon a constant 
background and have a width of approximately 25meV. 
As oxidation reactions are known to occur at a specific 
energy and at a length scale on the order of a , these 
structures represent a quantum dot structure. Advanced 
diagnostics to observe these quantum dots and the dy- 
namics of Redox reactions are topics of our current re- 
search and will be discussed elsewhere. At present we 
conclude that we have verified that Gauss' law applied 
to metal surfaces is an appropriate theoretical descrip- 
tion of surface charge which corresponds to the quantum 
mechanical model of a two-dimensional free electron gas. 

In Fig. 8.1b-8.5b we have plotted the observed density 
of states of charge under the condition that the surface 
of the metal plate is illuminated by UV light. The inten- 
sity of the UV source at the surface of the plates yields 
24mW/cm 2 in the case of Zn, Cu and Al and 96mW/cm 2 
for C and n-Si. The spectrum of the UV source that has 
been used is given in Fig. 6. As a result of the emission 
of electrons by the photoelectric effect, we observe a de- 
cline in the occupation level of the surface states. This 
decline corresponds to an increase in the occupation level 
of continuum states. 

The density of states of space charge is given by [18-21] 



dn 3d (E)/dE = V 2fc „ ^E-E . (12) 

The characteristic energy level E Q relates to the thresh- 
old energy level at which electrons are considered being 
in the continuum. 

We have fitted the data sets in Fig. 8.1b-8.5b (dashed 
lines) to the energy distribution function of the injected 
charge dQ/dE given by: 

dQ(E) _ A E<E Q 

dE \^A-^£n^E-^W E>E 

where A respectively Q define an effective surface area 
and volume. We have used E Q as a fitting parameter to 
obtain an estimate for the onset energy of space charge. 
We have fitted the line shape at the energy domain as 
indicated in Fig. 8.1b-8.5b. At higher energies (> 20eU) 
the temperature of the electrons in the continuum in- 
creases as a result of the ionisation of molecular nitrogen 
(data not shown). This results in an overestimate in the 
density of states in the high energy limit as predicted by 
Eq. 2 and 13. By using the fitting procedure as discussed 
we obtain a proper, unbiased fit of all our data to a sin- 
gle line shape. The threshold energies that have been 
determined are presented in Table III. 
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FIG. 8: The charge density on metal plates as a function of the Fermi level, a) In absence of UV light the density of states is a 
constant, b) In presence of UV light the surface charge is emitted and the density of states is a combination of surface charge 
and space charge. The dashed lines are the fitted shapes according to Eq. 13. It represents the transition of a surface bound 
electron gas into the continuum. 
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TABLE III: Observed threshold energies for the photoelectric 
effect determined by the method that has been described in 
this work compared to the Hartree-Fock estimate for exchange 
energy [18]. 



Element (a /r s 


) Observed threshold 


Hartree-Fock 


energy (eV) 


Exchange energy (eV) 


C 0.01 


0.2±0.4 


0.1 


n-Si 0.32 


3.5±0.2 


4.0 


Cu 0.37 


3.8±0.2 


4.6 


Zn 0.43 


4.5±0.2 


5.3 


Al 0.48 


5.7±0.2 


6.0 



V. DISCUSSION 

The observation of space charge, and more specific, the 
observation of an electron energy distribution that scales 
as VE, has been reported already in 1911 by CD. Child 
[28]. Child studied the emission current density from 
hot filaments in vacuum tubes and concluded that the 
current is limited by the presence of space charge. The 
phenomcnological description of this effect has led to the 
Langmuir- Child Law (or the three-and-a-half-power law). 
Upon differentiation of this law we obtain a density of 
states that is proportional to \[E. No reports have been 
found in literature where an offset energy is observed. 

More recent work deals with cold emission processes 
of electrons held in nano structures [1-3]. The observed 
energy distribution found in our work is certainly dif- 
ferent from electrons that result from high field field 
emissions [1-3,18,19,24]. In the latter case electrons are 
emitted at high electrical field strengths ranging from 10 
to 50kV/cm. The energy distribution follows a Fowler- 
Nordheim law as a result of quantum tunneling through 
the cathode-anode barrier [1-3]. Emissions that have 
been observed in our work occur in the energy range of 
to 5V, which is nearly two orders of magnitude less than 
reported in [1-3]. To ascertain that we actually deal with 
a different emission process, we studied the onset of field 
emission in our set-up. In order to observe electron cur- 
rents resulting from this process we had to increase the 
electric field strength to a level in excess of 30 kV/cm. 
This was realised by reducing the distance between the 
metal plate and the grid to 5mm and switching to an 
operating voltage of 15kV (data not shown). 

We studied the onset energy for the photoelectric effect 
for a range of metals: Zn, Cu, Al, C and n-Si. In all cases 
we observed a positive threshold energy. Zn has a work 
function of 4.2 eV, a value that is accurate within 5% 
[19]. We studied the onset energy for Zn by varying the 
UV photon energies from 6.7eV to 4.5eV. We found that 
the onset energy for emission at +4.5 eV does not change 
within the experimental error of ±0.2eV. By analysis of 
the onset energy found for Zn, Cu and Al (Table III) we 
can conclude that emission occurs in general when the 
surface is charged and the electrons reach an energy level 
associated with the work function of the metal. The spe- 



cific behaviour of the Fermi-Dirac distribution, which has 
been discussed in Fig. 3, accounts for a sharp increase 
in the amount of space charge when the Fermi level ex- 
ceeds the continuum level. Thus, by statistical analysis 
alone, an accurate quantitative prediction has been made 
for the work function which is determined by the onset 
energy for space charge. 

The observed emissions cannot be explained by the 
conventional interpretation of the photoelectric effect. In 
the latter case, the threshold energy for the emission of 
photoelectrons for Zn is expected at <p — hv — —2AeV 
when using 180nm (6.7eV) photons [26-27]. As a result, 
these observations let us no other choice than to con- 
clude that the observed emissions are of different origin. 
This conclusion may be implausible from a historic point 
of view. However, the Einstein-Millikan threshold rela- 
tion hv — <f) accounts for the minimum photon energy 
based on the measurement of the kinetic energy of elec- 
trons upon arrival. This is different in our experiments, 
where we control the Fermi level of an electron ensem- 
ble and observe when emission occurs. With certainty, 
electrons that satisfy the Einstein-Millikan relation arc 
also released in our experiments. However, these ejected 
electrons have energy of a few eV above the Fermi level. 
They therefore loose their kinetic energy in collisions with 
gas molecules that are present in the ambient environ- 
ment leading to a recapture by the surface. As a result 
these electrons do not contribute to the process of elec- 
tron emission. It is by the construction of the experiment 
that we study the fate of the electron system as a whole. 
This is in contrast with the stopping voltage method [26] 
which is designed to determine the photon energy by ki- 
netic analysis of the photo electron. It should not come 
as a surprise that the physical processes involved are of 
different origin. At this stage we conclude that our re- 
sults do not contradict the Einstein-Millikan relation. 

The Fermi level in this work is controlled by electro- 
static charging. As a consequence, the emission of elec- 
trons by the photoelectric effect is not limited by en- 
ergy requirements. Indeed, we have observed an onset 
in emission at the point where the Fermi level coincides 
with the continuum level. We conclude that nearly free 
electrons are involved in the emission process. This im- 
plicates that the ejection is the result of a binary process, 
involving one photon and one electron. We consider the 
simplest process: an elastic collision of a photon and an 
electron (Compton scattering). The process is schemat- 
ically depicted in Fig. 9. As energy is conserved in the 
scattering process, the loss in the photon energy is equal 
to the kinetic energy Ek that is gained by the electron. 
In a head-on collision (9 = tt) this yields Ek = h 2 k 2 /2m, 
where fe is the wave vector of the incident photon. For 
visible light, the energy that is exchanged in the process 
is on the order of 10/ieV. This is insignificant when com- 
pared to the amount of thermal energy that is available at 
300K, which yields kT=26 meV. The amount of momen- 
tum that is transferred to the electron is considerable. It 
equals twice the photon momentum i.e. p — 2hk. The 
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i, =\ + — (l-cos<?) 
mc 

FIG. 9: Elastic scattering of electrons by photons (Compton 
scattering). Left) A photon approaches an electron Right) 
The electron is scattered and has gained a momentum in the 
process. The energy and momentum exchanged is determined 
by the Compton shift Ai — Ao. 



average amount of momentum that is exchanged by the 
Compton process equals p — hk — 2irh/X. We argue that 
this corresponds to the minimum amount of momentum 
that is needed to escape from a surface. 

For an electron to escape a surface momentum is 
needed. This is the consequence of the broken transla- 
tion symmetry near the metal-air interface, in compliance 
with Noethers' theorem on conservation laws and sym- 
metry [21]. In order to satisfy the conservation of total 
momentum, a source is required that delivers an electron 
a sufficient amount of momentum at emission. We shall 
estimate the minimum amount that is needed by evalua- 
tion of the kinematic aspects of the ejection process. 

Consider the momentum eigen states p = Hn/a(i,j,l) 
of an isotropic free electron gas, given by a set of quantum 
numbers {i,j,l) that corresponds to the energy states 
Eiji = h 2 k 2 /2m [i 2 + j 2 + I 2 ] where k stands for the 
magnitude of the wave number associated to the elec- 
tron. This wave number is determined by the quanti- 
sation length a through k = ■n/a. For surface charge 
i.e. an electron gas with two degrees of freedom, l~ 
and only (i, j) states are occupied. For an electron in 
the continuum no restrictions have to be made and all 
states have to be considered. This means that in 
a transition from the surface into the continuum a mini- 
mum amount of momentum is required perpendicular to 
the interface of magnitude p = hk = h/a. The quantisa- 
tion length a is specified by the local symmetry of space 
near the surface that is defined by the surface structure 
of the lattice at the metal-air interface. Structure is de- 
fined either by geometrical constraints or by the presence 
of impurities. In our experiments the surface structure 
of the metals and the graphite is specified by the micro- 
scopic surface roughness which lies in the range of 0.1-1 
/im. From the dopant concentration in our n-Si sam- 
ple (~ 10 15 cm~ 3 ) we deduce that the average distance 
between impurities yields ~ 0.1 /j,m. Thus, the quantisa- 
tion length of the surface structures in our experiments 
is in the order of the wavelength of UV light. We believe 
that the observed threshold wavelength of the light Xthr, 
corresponds to the threshold of electron momentum n/a 
that is needed to remove an electron from the surface: 
\ t hr — 2a. This is in agreement with our observations 



that the threshold wavelength does not depend on the 
Fermi energy. For this, we have charged Zn plates to 
-100V, which is far beyond the continuum level and illu- 
minated the surface with intense sources of light (UV-B, 
UV-A and visible light), all of well defined wavelength. 
We have probed for the onset of emission and found that 
the threshold wavelength to observe electron emission lies 
below 300nm. We therefore conclude that under the ex- 
perimental conditions presented in this work the ejection 
of electrons is triggered by Compton scattering. 

For carbon we have found a value for the work func- 
tion that yields 0.2±0.4 eV (see Fig. 8.4b and Table 
III). This value is consistent with our assumption that 
exchange energy is the major contribution to the work 
function. Graphite is a semi-metal with a free electron 
density several orders of magnitude lower that in met- 
als [18,19, Table I]. Exchange mechanisms are therefore 
absent, and as a result, the work function is non-existent. 

Silicon differs from the other elements under consider- 
ation since the sample that was used is crystalline. Due 
to the hardness of the material its surface could not pre- 
pared to obtain a polycrystalline sample. No attempts 
where made to remove any residual oxides from the sur- 
face. We have found structures in the surface density of 
states in the region of 0-5eV (Fig. 8.5a and 8.5b) which 
we attribute to band structure effects in the bulk and 
surface decontamination. Despite the fact that the sili- 
con data in Fig. 8.5b deviates from the fitted model in 
the region of 0-5eV we have included the estimate for the 
threshold energy in Table III. This is motivated by the 
fact that a detailed, microscopic understanding of these 
effects is not a prerequisite for a correct description of 
the macroscopic parameters that describes a many par- 
ticle system [20]. In this work we use statistical physics 
to describe the macroscopic parameters surface charge, 
Fermi energy and exchange energy. The use of statistical 
physics guarantees that the qualitative and the quanti- 
tative behaviour of a system is well described, at least 
in first order, independent of microscopic details. As our 
statistical model discribed by Eq. 13 fits to the overall 
shape of the Si data and since the threshold is a rea- 
sonable estimate of the work function assuming a bulk 
electron density for phosphorus doped n-Si, we have in- 
cluded Si in our work. 

Based on all of the results presented in Table III we 
conclude that the work function of a metal is linearly 
dependent on the electron density l/r s according to 
W = (11.0 ± 0A)a o /r s . This is consistent with the scat- 
tered work function data found in literature, presented in 
Fig. 1, which yields W = (9.9)a /r s on average. As the 
Hartree-Fock theory predicts that the exchange energy 
is a linear function in electron density we conclude that 
the bonding of electrons to a surface is dominated by the 
process of electron-electron interaction. The remaining 
discrepancy between the experimental values found and 
the Hartree-Fock estimate is attributed to screening. The 
effect of the screening of positive ions by strongly bound 
electrons is neglected by the model, but tends to reduce 
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the exchange energy by 10-20% [18]. 



By taking into account the sharp onset in the energy 
distribution that has been observed in Fig. 8.1b-8.5b, 
we suggest that the photoelectric effect is the result of 
phase transition from a condensed phase (surface charge) 
to a vapour (space charge) . This is motivated by the ob- 
servation that the transition can be crossed only in one 
direction (initially an excess amount of negative charge 
is needed). Secondly, we have found that the sharp on- 
set coincides with the energy needed to break electron- 
electron bonds. This occurs when the Fermi level of an 
electron ensemble held at a surface, has been risen by an 
amount that matches the work function. It results in a 
phenomenal increase in the amount of space charge by a 
factor of e qV ' kT — e 162 which yields 70 orders of mag- 
nitude. This abrupt change in the occupation number 
of electrons in the continuum states (i.e. evaporation) 
is similar to the behaviour of the distribution function 
required for bosons in the low temperature limit. The 
number of bosons in the ground state increases abruptly 
(i.e. condensation) after a thermal equilibrium is reached 
below the Bose-Einstein transition temperature [20,21]. 
In the case of a dense electron system thermal equilib- 
rium is only reached at temperatures that are in the or- 
der of the Fermi temperature. This associated temper- 
ature typically yields T=20,000 K for metals. At ther- 
modynamical equilibrium with black body radiation at 
these extreme temperatures, the average photon energy 
of E = hc/X = 2.3 10~ 4 T=4.6eV and the photon den- 
sity of n = 3.0765/A 3 w 10 20 m~ 3 [20,21] are certainly 
sufficient to account for emissions. This corresponds to 
the high temperature limit where thermo-ionic emissions 
(T~ 2000 K) are observed [18,19]. In the low tempera- 
ture limit black body radiation does no longer account for 
emissions. We found, however, that the presumed ther- 
mal equilibrium can be effectively mimicked by raising 
the Fermi level into the continuum and illumination of 
the surface by an auxiliary source of UV photons (Fermi- 
Dirac evaporation). 



Finally, we argue that the photoelectric effect is an in- 
duced phase transition i.e. a process that is triggered 
after applying a force by using an analogy with super- 
heated water. In the latter case, the temperature of the 
liquid phase exceeds the temperature of the vapor phase. 
Only after applying force to the water surface, the van- 
derWaals bonds between the water molecules are broken. 
This results in a sudden release of water vapor leading to 
the immediate cooling of the liquid. In this way the ther- 
modynamical equilibrium between the liquid and vapor 
phase is restored. This is similar to surface charges held 
by a conductor at a negative potential; its excess charge 
is abruptly released after introducing photons which me- 
diate forces on electrons. 



VI. CONCLUSIONS 

We have presented and validated a method to measure 
the energy distribution of dense electron ensembles at 
ambient conditions. The method relies on the direct mea- 
surement of the density of states at the Fermi level. By 
applying Bohrs' correspondence principle we have found 
that surface charge, as treated by Gauss' law in classi- 
cal electrostatic theory, is the macroscopic limit of a free 
quantum electron gas with two degrees of freedom. The 
Fermi level charge held at the surface of a metal plate is 
determined by its electrical potential through E = —qV. 
When the Fermi level matches the work function, a sharp 
onset in the density of states is observed for Zn, Cu, Al, 
n-Si and C (as graphite). The structures correspond to 
the onset of space charge i.e. electron emission at respec- 
tively 4.5±0.2, 3.8±0.2, 5.7±0.2, 3.5±0.2 and 0.2±0.4 
eV. The characteristic onset energy equals the work func- 
tion. Remarkably, the work function of Carbon (a low 
electron density material) is vanishing. After compari- 
son of the experimental values for the work function to 
the Hartree-Fock model of a free electron gas we conclude 
that exchange energy is the major factor that contributes 
to the work function of metals. 

Electron emissions were triggered in our work by the 
photoelectric effect. However, the observed emissions in 
this work do not obey the Einstein-Millikan threshold 
relation hv = <j> of the photoelectric effect. By discrimi- 
nation between requirements for energy and momentum, 
we find that the observed emissions in this work are trig- 
gered by photon momentum. Therefore we conclude that 
the photoelectric effect is the result of a phase transi- 
tion in an electron gas that connects a surface state with 
two degrees of freedom to a continuum state with three 
degrees of freedom. The assignment of a phase transi- 
tion to the observed effect has been justified by analyz- 
ing the behaviour of the Fermi-Dirac energy distribution. 
When the Fermi level is risen above the continuum level, 
the amount of space charge increases abruptly. The in- 
crease in the number of continuum electrons is phenome- 
nal and exceeds 70 orders of magnitude. The experimen- 
tal threshold energies and the line shapes of the transition 
are in good agreement with theoretical predictions. 
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